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Abstract 

We study the magnetic excitation spectra of L-edge resonant inelastic x-ray scattering (RIXS) 
from the spin singlet ground state in one-dimensional undoped cuprates. Analyzing the transition 
amplitudes of the magnetic excitations in the second-order dipole allowed process, we find that the 
magnetic excitations are brought about not only on the core-hole site but also on the neighbor- 
ing sites. The RIXS spectra are expressed by the one-spin correlation function in the scattering 
channel with changing polarization, and the two-spin correlation function in the scattering channel 
without changing polarization. The latter could not be brought about within the so-called UCL 
approximation. We calculate these correlation functions on a finite-size ring. An application to the 
possible RIXS spectra in Sr2Cu03 demonstrates that the contribution of the two-spin correlation 
function could be larger than that of the one-spin correlation function in the a polarization for the 
momentum transfer around the zone center. 

PACS numbers: 78.70.Ck, 72.10.Di, 78.20.Bh, 74.72. Cj 
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I. INTRODUCTION 



Resonant inelastic x-ray scattering (RIXS) has attracted much interest as a useful tool to 
investigate excited states in solids. Both the i^-edge and the L-edge resonances are available 
in transition-metal compounds. To observe the momentum dependence of the spectra, it is 
suitable to use the K-edge resonance that the ls-core electron is prompted to an empty 4p 
state by absorbing photon and then the photo-excited 4p electron is recombined with the 
core hole by emitting photon. It is known that the intensities for energy loss u above several 
eV's arise from charge excitations to screen the core-hole potential.- - — A useful formalism 
has been developed^— on the basis of the Keldysh Green function, 15 in which the RIXS 
spectra are described in terms of the 3<i-density-density correlation function. The spectra 
have been calculated for undoped cuprates,— ~— NiO,— and LaMn03, 17 in good agreement 
with the experiments. In addition to such spectra, the intensities for u below 1 eV have been 
observed from the Lf-edge RIXS in L^CuO^.— ^ Their origin is attributed to the magnetic 
excitations caused by the modification of the exchange coupling under the presence of the 
core-hole potential.— >2i The above-mentioned formalism has been adapted to describe such 
magnetic excitations, having led to the expression that the RIXS spectra are proportional to 
the two-spin correlation function, which has been evaluated by the 1 / S'-expansion method.— 

Recently, the L-edge RIXS experiments with high resolutions have been carried out in 
transition-metal compounds. The observed spectra for u around several hundreds meV's 
can be understood as coming from magnetic excitations.— ~— The L-edge resonance in the 
undoped cuprates is described by the process that the 2p-core electron is prompted to the 
empty x 2 — y 2 orbital by absorbing photon, and then an occupied 3d electron combines with 
the core hole by emitting photon. If the 3d orbital in the photo-emitting process is different 
from the one in the photo-absorbing process, the excitations within the 3d states are brought 
about ("d-d" excitation). 26 Even if the 3d orbitals in the photo-absorbing and photo-emitting 
processes are the same, purely magnetic excitations could be generated when the direction of 
the staggered moment deviates from the z direction.— The spectra observed as a function of 
energy loss exhibit a systematic variation with changing momentum transfer.— That is, the 
peak position moves according to the spin-wave dispersion curve, while their shapes exhibit 
structures indicative of two- or three-magnon excitations. 

In our previous paper,— we have analyzed the process generating the final states in 
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the second-order dipole allowed process, and have clarified how the spin-flip excitations are 
brought about on the core-hole site. We have also found that the spin-conserving excitations 
are brought about around the core-hole site due to the strong perturbation working in the in- 
termediate state. Note that, within the conventional approach called UCL approximation,— 
no spin-conserving excitation is possible, and the spin-flip excitation is generated just on 
the core-hole site. Since no core hole exists in the final state, these excitations are freely 
moving in crystal. We have treated the quantum fluctuation in the final state within the 
1/S expansion method, which is known to work well in the two-dimensional 5 = 1/2 Heisen- 
berg antiferromagnet.— On the basis of these results, we have analyzed the RIXS spectra in 
Sr 2 Cu02Cl2. The result shows that substantial intensities are generated in the high energy 
side of the one-magnon peak, which are originated from two- and three-magnon excitations, 
in good agreement with the experiment.— 

The above analysis has been carried out in the antiferromagnetic ordered state. Although 
the spin-flip excitations are likely to extend on neighboring sites due to the second-order 
process, the systematic analysis of such effects is quite difficult there, because the spherical 
symmetry is broken in the spin space. On the other hand, in the one- dimensional Heisenberg 
antiferromagnet, since the ground state is the spin singlet, the spherical symmetry remains 
intact in the spin space. Such a system is quite suitable to study how the spin excitations are 
generated around the core-hole site. The purpose of this paper is to clarify the nature of the 
magnetic excitations in RIXS by studying the spectra on the one-dimensional system. We 
analyze the second-order process on a microscopic model, with distinguishing the spin-flip 
and spin-conserving excitations. As a result, the RIXS spectra are expressed by the one-spin 
and two-spin correlation functions in the scattering channels with and without changing po- 
larization, respectively. Note that the former includes the effect of the magnetic excitations 
generated on the neighboring sites. We could evaluate convincingly the amplitudes leading 
to these correlation functions on a finite-size ring. The amplitudes of the excitations on 
the neighboring sites are found to increase with decreasing values of the life-time broaden- 
ing width of the core hole. The present analysis contrasts with the previous analysis by 
Haverkort,— where the form of the magnetic excitations is inferred from Hannon's formula^ 
of the resonant elastic scattering. The analysis as well as the UCL approximation applied 
to the one-dimensional system^ could not describe the extended nature of the excitations. 
Finally we analyze possible RIXS spectra in Sr 2 Cu03, demonstrating that the contribution 



3 



of the two-spin correlation function becomes larger than that of the one-spin correlation 
function in the a polarization, and could be distinguished by comparing the spectral shape 
as a function of energy loss between in the a and it polarizations. 

The present paper is organized as follows. In Sec. [TTJ we describe the Hamiltonian 
responsible for magnetic excitations and transition-matrix elements relevant to the L 2j 3 edges 
in cuprate compounds. We analyze the second-order process giving rise to the magnetic 
excitations through the intermediate state, and express the RIXS spectra in terms of spin- 
correlation functions. The amplitudes leading to the spin-correlation functions are evaluated 
on a finite-size ring. In Sec. IIII[ we numerically calculate the spectra on a finite-size chain. 
Section [TV] is devoted to the concluding remarks. In Appendix, the L 2>3 -edge absorption 
coefficient is briefly discussed. 

II. FORMULATION OF RIXS SPECTRA AT THE L 23 -EDGE 
A. Second-order dipole allowed process 

Aiming at the application to the one- dimensional cuprates such as Sr 2 Cu03 and SrCu0 2 , 
we consider each Cu atom has one hole in the x 2 — y 2 orbital at the half-filling, where the x 
and y axes are defined along the Cu-0 bonds and the z along the crystal c axis. The low- 
energy spin excitations are described by the one- dimensional antiferromagnetic Heisenberg 
Hamiltonian, 



with 5=1/2 and J > 0. The ground state is known to be the spin singlet due to the 
quantum fluctuation. Unlike the antiferromagnetic ordered state, the spin singlet state 
possesses no special direction in the spin space. To manifest the rotational invariance in the 
spin space, we define the spin coordinate frame of x', y', z' axes by rotating the crystal- fixed 
coordinate frame of a, b, c axes with the Euler angles a, (3, and 7. 

In the electric dipole {El) transition, a 2p-core electron is excited to the 3d states at the 
transition-metal L 2j 3-edge. The 2p states are characterized by the total angular momentum 
j = 3/2 and 1/2 due to the strong spin-orbit interaction. The eigenstates with j = 3/2 



may be expressed as |& |), y/TjZ^ I) + VV^l^o t), v^l^o I) + a/V^-i t), 10-1 I), 



for m = 3/2, 1/2, —1/2, —3/2, respectively, and those with j = 1/2 may be expressed 




(2.1) 
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as -^2/310! |) + v/VS^o t>, "V^l^o I) + v^l^-i t>, for m = 1/2 and -1/2, 
respectively, where m represents the magnetic quantum number. The orbitals <f>x, <po, and 0_i 
have the same angular dependence as the spherical harmonics Y\i, Y w and Yi_i, respectively. 
The coordinate frame for both orbital and spin of the core hole is fixed to the crystal, which 
is different from that for the spins of the 3d states, that is, and 4- are associated with the 
direction of the crystal c axis. The corresponding interaction between photon and electron 
at site i may be described as 

1 



#int = w V — == V D»(jm, cr)h\ m c^d la e ici r ' + H.c. 



(2.2) 



where c qM stands for the annihilation operator of photon with momentum q and polarization 
li. The hj m represents the creation operator of the 2p hole with jm, and di a denotes the 
annihilation operator of the 3d hole with the x 2 — y 2 orbital and spin a. The w is a constant 
proportional to J °° r 3 R 3 d(r)R2 P (r)dr where i?3d(r) and R2 P (r) are the radial wave-functions 
for the 3d and 2p states of Cu atom. The D M (jm, a) describes the dependence on the 
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core-hole state and 3d spin, which is shown in Table I of Ref. 

With this interaction, the RIXS spectra may be expressed by the second-order dipole 
allowed process: 

2 

(<$> f ,\H int \n)(n\H int \<$>i) 



W(q f ocf,q, 
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x 5{E g + 



E 



Eg + Ui 

Eft - u f ), 



Em 



(2.3) 



with q { = (cu,Wi), q f = (q/,w/), |$<) = c^g), \® r ) = c qfCtf \f), where \g) and \ f ) repre- 
sent the ground state and excited states of the matter with energy E g and Ef, respectively. 



B. Magnetic excitations around the core-hole site 

Assuming that the core hole is created at the origin in the intermediate state, we write 
the ground state \g) of H mag as 

l<7>HtM) + UM>, (2-4) 

where | f) and | \) represent the spin states at the origin, and and \i)jq) are constructed 
by the bases of the rest of spins. 
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Then, just after the El transition, the wave function may be written as 



£ D^(jm,a)\r ) 



\jm), 



(2.5) 



where \jm) represents the core hole state. Note that the spin degree of freedom is lost at 
the origin in the intermediate state. 

Let H' be the Hamiltonian in the intermediate state. The states [iffy and \i[)q) are not 
the eigenstates of H' . Introducing the normalized eigenstate |0^)'s of H' with eigenvalue e' 
we rewrite the second-order optical process as 

1 
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*r - el 



(2.6) 



(2.7) 



+ e g ~ e co: 

where e g represents the ground state energy of H ma , g . The e core denotes the energy required 
to create a core hole in the state \jm) and the 3(i 10 -configuration. The T stands for the 
life-time broadening width of the core hole; T ~ 0.3 eV at the Cu L 3 edge. We write the 
polarization-dependent factor by 



m 



(2.8) 
(2.9) 



Table [J shows P^ and P^ for on and aj along the x, y, and z axes, where sgn(a) = 1 
for u —"\ and —1 for a =\,. 



Scattering channel with changing polarization 

We analyze the scattering channel that changes the polarization during the process at 
the L 3 -edge. Let an and «/ be along y and x axes, respectively. We have the spin-conserving 
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TABLE I: PP \j ;atf,oti ) and P a '{j \af,ai) where upper and lower signs correspond to a =f and 
X, respectively. 
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term coming from Pi°^ and the spin-flipping term coming from Pp. The spin-conserving 
term is given by 



E I n ^ . i f - F ^ ' ^ int ' 5 ^ 



&i + Eg — E n 

« n (0) i t) E i^(<k^iv4) + n (0) i E i^>*«)toM>> 



7J JJ 

Z 



cx ( -- ) cos/3 { | t) E 1^)^(0(^14) " I I) E 1^)^(0(^1^) \ • (2-10) 



We consider the states Sq \g) and (S(' + S 2 L 1 )\g) onto which Eq. (12.1 Oft is projected. Although 
they are orthogonal to \g), they are not orthogonal to each other nor normalized. Let 
and |^ 2 ) be Sq \g) and (S" +Sti)\g), respectively. Then, the overlap matrix (p z ')i,j = (V'ilV'i) 
is given by 

7 2a 



■1 



2a i + 26 



(2-11) 



where a and 6 are static correlation functions between the nearest and the second-nearest 
pairs, respectively, that is, a = (S* Sq) = (S^Sq), b = (Sf S^}. Using the inverse of p z >, 
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we obtain 



5> 



mt\n) 



+ Eg — E n 



(n\H iQt \g) 



oc 



^) cos/3 {./ffa)^' + / 2 {1) (^)(5f + S*)} |«?>, (2.12) 



where 



with 



(2.13) 



ra=l,2 



(^)i = (^lEl^K)^), 

(Q«) 2 = 4(^ Vf E 1^)^(0(^1^) ■ 
Therefore, the final state is expressed as 

cos/3 L/fVOSS' + f?(^)(S( + 5^)1 |«?) 



(2.14) 
(2.15) 



i 

15 



-— ] a f± x Oi ± • 



/ 1 (1) (a; l )So|| + / 2 w (a; l )(S 1 ||+S_ 1 ||) \g) (2.16) 



c(l) 



where Sjy stands for the z' component of Sj. The etij_ and otf± represent the polarization 
vectors projected onto the x-y plane. Therefore etf± x oci± is always parallel to the z axis. 
The spin-flip term may be expressed as 

1 



5> 



int w 



+ E„ — E n 



(n\H int \g) 



oc (±) sin/3 je-*| |) E + ^\ t> E l^>««X^I*| ■ ( 2 - 17 ) 

We project this state onto a pair of states S^lg) and (S^ + Sl-JIg), and another pair of 
states S^lg) and (S^ + SZi)\g), where Sf = Sf ± iS? . The states belonging to different 
pairs are orthogonal. Defining the overlap matrices p + and p_ for each pair in a similar way 
used for p z >, we immediately notice that 



P+ = P- = 2pz 



(2.18) 
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The second term of Eq. (12.171) is projected onto the first pair of the states, leading to 



v | ^) sin^ {g 1 (co, l )S^ + g 2 (co, l )(S+ + S^)} \g), (2.19) 



where 



with 



9m(Ui) = (P+UniQ+U ( 2 - 20 ) 

n=l,2 



(Qi ] )i = (4\ e i^R^wt) = (Q^h (2.2i) 

v 

(Q«) 2 = 2<V>J|3r £ |0,>i2«)(^|V*> = (Q?*)* (2.22) 

Hence we have gi(uJi) = (1/2) /i (wj) and g 2 (uJi) = (1/2) / 2 (^i)- The similar relations are 
obtained for the first term of Eq. (12. 17ft by using the second pair of the states. Therefore 
Eq. (I2.17P is expressed as 



oc 



sin/3{cos7[/i (1) (wi)S? + fP(<»i)(Sf + S^)) 

= «A x «*± • [A (1) (^)So± + / 2 (1) N)(Si± + S_ij_)]|flr>, (2.23) 

where Sj± stands for the component perpendicular to the z' axis of S^. Combining this 
result with Eq, (I2.16p . we finally obtain 

V H int \n) — — (n | H int \g) 



-- 1 a f xtxi- [/r^^So + / 2 w (wO(Si + S_i)]b). (2.24) 

Excitations on further neighboring sites are included by adding a term /^(o^)^ + S_ 2 )|5') 
with f^\uji) similarly determined coefficient. 

Collecting up the amplitudes from all sites, we obtain the expression of the RIXS spectra 
for the polarizations a^f) = (af (/) , a^ (/) ), 

4 / J \ 2 

^(g/,^,^;^,^,^) = ^— J (o^aV - a y f a^) 2 Y {t \uji-q,uj), (2.25) 



where 



F«(^; g ,a;) = J (Z^(co l ;q,t)Z^\u; i ;q,0))e iuit dt, (2.26) 
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with 

ZW( Wi ; q) = YllfP^S; + / 2 (1) (^)(^ + i + + / 3 (1 V*)(^ + 2 + ^- 2 )]e"^ 

= [fi 1 \u l ) + 2^ 1 \u t )cosq + 2jf\u l )cos(2q)}Y,S]e-^. (2.27) 

The presence of f!^\ui) and /^(wj) modifies the g-dependence of the spectra, but the 
spectral shape as a function of u is expressed by the conventional correlation function. 



2. Scattering channel without changing polarization 

In this scattering channel, only the spin-conserving excitations are brought about through 
the diagonal components of Pj?\ For each component, we have 

Yl Hint i n ) , , i J" — w ( n i Hint i 9 ) 

LOi + hi „ — hi n 

n ^ 

x (YE) {' T) E I^^O^I^J) + I E I6»>*«MI^>} • (2-28) 

Since the polarization dependence behaves as otf±-cti±, the corresponding spin excitations 
may be expressed by operating even number of spin operators on the ground state. Therefore 
it is reasonable to assume the excited states as (Si + S_i) ■ S \g), (S 2 + S_ 2 ) • Sols'), an d 
Si • S-ilp). Since these states are not orthogonal to \g), we define the overlap matrix p 
by including \g) to be projected in addition to the above three states. The procedure of 
projection is the same as before. Here we simply write down the result; 

Yl Hmt ' ^ T J" F~ ^ I Hint 1 9 ) 

^-^ LUi + hi q — hi n 

oc ^a f± ■ Oi^fPiui) + / 2 {2) (^)(Si + S_i) • So 

+ /f(o; l )(S 2 + S_ 2 ) • So + /f(^)Si • S_i}|^). (2.29) 

Collecting up the contribution from all sites, we obtain the RIXS spectra for the polar- 
ization vectors cti and ctf, 

W(qf,Uf,a f ]q i ,u i ,a i ) = -^-\^—J (a x f a* + a)(x(f F (2) (^;g,w), (2.30) 



where 



Ui; q, t)Z^ 2 \ui\ q, 0))e ZOJt dt, (2.31) 
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with 

Z^\u i]q ) = J2{fP (w i )(S i+1 + S i _ 1 ).S i + / s (2) ( Wi )(S i+2 + S i _ a ).S i 
j 

+ /^(wOSj+i-S^iJe^. (2.32) 

In the far-off-resonance condition that |cjj + e 9 — e cor c| 3> I\ e' and in the UCL approx- 
imation that T ^> \uji + t g — e core |,e', we could factor out R(e') in Eq. (12.28 p . In such 
circumstances, using the relation Y2 V I^K^I = 1, we notice that Eq. (I2.28P is proportional 
to \g), that is, no excitations are generated. 

C. Evaluation of coefficients 

For evaluating /^(w^'s and (u;j)'s, we consider a system consisting of 12 spins of 
5 = 1/2 with periodic boundary conditions for the initial and final states, as shown in Fig. [TJ 
Since the relevant magnetic excitations are restricted around the core-hole site, we expect 
that a system having rather small size works well. Representing if mag by a matrix of 924 x 924 
dimensions in the subspace of the z' component of the total spin S^ ot = 0, we diagonalize the 
Hamiltonian matrix. We obtain the ground state energy as e g / {NJ) = —0.448, which should 
be compared with the exact value —0.443.— On the other hand, the intermediate state is 
expressed by using 11 spins, since the spin degree of freedom is lost at the core-hole site. 
Therefore H' may be represented by a matrix with 462 x 462 dimensions in the subspace of 

Sit = ±1/2- 

As briefly discussed in Appendix and shown there in Fig. |6l the absorption coefficient for 
the transition 2p — > 3d x 2_ y 2 in cuprates has a single-peak structure; the peak is located at 
Ui = cu® = e cor e + 0.8J. Note that the incident photon energy uji is usually tuned to give 
the maximum absorption coefficient in RIXS experiments. For Sr 2 Cu03 and SrCu02, J is 
rather large (200 — 250 meV), and V ~ 0.3 eV, and thereby we have Y / J ~ 1.2. 

Using the eigenvalues and eigenfunctions on a ring of 12 spins, we calculate the coefficients 
for Ui = Table HT1 lists the calculated values. For ffp's, If^l and \ f^ \ are rather smaller 
than \fi \ with Y/J= 1.2, while they become larger with Y/J — 0.5, indicating that the 
effect of magnetic excitations on neighboring sites increases with decreasing value of Y. As 

(2) (2) 

regards fh s, |/ 2 I overwhelms other absolute values even when Y/J = 0.5. This suggests 
that the disturbance is nearly limited within the nearest neighbor sites. 
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FIG. 1: A ring of 12 spins used to evaluate (wj)'s and fi, (wj)'s. The spin at site is annihilated 
in the intermediate state. 



TABLE II: Coefficients fji ViO's and /r^j's in units of 1/J. 



r/j 


A (1) (^°) 


/ 2 (1) K°) 


/ 3 (1) K°) 


1.2 
0.5 


(0.182,-1.752) 
(0.649, -4.395) 


(0.282,-0.188) 
(0.723, -0.979) 


(0.165,-0.106) 
(0.426, -0.567) 




r/j 


/i 2) K°) 


/i 2) K°) 


/i 2) K°) 


1.2 

0.5 


(0.038, -0.365) 
(-0.431,-0.543) 


(-0.007,0.009) 
(0.010,0.022) 


(0.039,-0.067) 
(-0.074, -0.142) 



III. RIXS SPECTRA 



Since the magnetic excitations could propagate in the crystal in the final state because 
of the absence of core hole, a small-size ring would not work well. In the following, we 
calculate the correlation functions on a ring of 16 spins from Eqs. (I2.26f) and ( 12.311) . Figure 
[2] shows (u;?; q, lo) and Y^(u^; q, cu) numerically calculated with Y/J = 1.2. The blue 
lines represent the des Cloizeaux-Pearson curve,— u = J(tt/2) sing, which is the lowest 
boundary of the excitation energy. These functions have already been obtained from the 
calculation on a finite-size system 3 - 2 - and also from the Bethe Ansatz solution. 35 The spectral 
shape of Y^ 2 \uj^] q, uS) as a function of u seems to have more weights at higher u than that 
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FIG. 2: (Color online) Correlation functions calculated on a ring of 16 spins, as a function of energy 
loss oj for various q- values; (a)Y^ q, uj)/N and (b)Y^ (a;- 1 ; q, to) /N, with oj® the photon energy 
giving rise to the peak in the absorption spectra, and T/J = 1.2. The blue lines represent the des 
Cloizeaux-Pearson curve, which is the lowest boundary of the excitation energy. 

Figure [3] shows the integrated intensities denned by 

I (1) (^;q) = jY^;q,co)^, (3.1) 

I {2) (^ q ) = Jy*\u,Uu)t£- (3-2) 

The 7^(a;°; q) vanishes with q — > 0, increases with increasing values of q, and remains finite 
with q — y tt. The presence of f% (oJi) and /| makes the g-dependence deviate from 
that of the dynamical structure factor. The deviation becomes conspicuous with V / J = 0.5, 
because of the increase of l/jj an< ^ Ifi The I^{uj®\q) is found one order of 

magnitude smaller than /W(u;°; q) around the zone center. 

Bearing in mind a possible application to Sr 2 Cu03, we demonstrate the importance of 
Y^iojf] q,co) by calculating the RIXS spectra on the L 3 edge under a typical scattering ge- 
ometry shown in Fig. HI the same geometry as used in the experiment in S^CuC^C^;— the 
angle between the incident and the scattered x-ray is kept 130 degrees, and the scattering 
plane includes the b(x) and c(z) axes. The polarization vector of the incident photon is then 
expressed as a; = (0, —1,0) for the a polarization and oci = (Xi,Q,Xi) f° r the 71 polariza- 
tion. Similarly, the polarization of the scattered photon is expressed as otf = (0, —1,0) for 
the a' polarization and olj = (xj,0,X/) for the tt' polarization. The polarization is usually 
separated with the incident photon, but not separated with the scattered photon in experi- 
ments. In such a situation, we may express the RIXS spectra depending on the polarization 
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ll/2 

q 



n/2 

q 



FIG. 3: (Color online) Frequency-integrated intensities of the correlation functions, (a) 
I^(uj°-q)/N and (b) 1^ (w° ; q) / 'N , calculated on a ring of 16 spins as a function of q. The 
lo® is the photon energy giving rise to the peak in the absorption spectra. Filled black and open 
red circles correspond to Y/J = 1.2 and 0.5, respectively. Curves are guides to the eye. 



of the incident photon as 



w 



x 



YW(u t] q,u)+{±) 2 YM(u t ;q,u) 

LUi] q,u) 



2 Y^q,u) + (^f) 2 Y^ 



[a 



7T 



pol. 
pol. 



(3.3) 

where q is regarded as the transferred momentum projected onto the b axis. The contri- 
bution of Y^iuji] q, to) relative to that of Y^\uji; q,oj) is enhanced by {2/x 7 }) 2 in the a 
polarization. The contribution of Y^(ui, q, to) in the tt polarization is reduced from that in 
the a polarization by a factor (x/X7) 2 - 

Figure [5] shows the RIXS spectrum as a function of energy loss u for q = 3tt/A(— 37r/4), 
where xl = 0.89(0.23), x} = 0.23(0.89) for u { ~ 930 eV and b = 3.49 A. We put J ~ 260 
meV and Y/J = 1.2. The calculated curves are convoluted by the Lorentzian function 
with the half-width-half-maximum of the possible resolution, 78 meV. The contribution of 
yd) 

(uf;q,uj) dominates the spectra in the n polarization at both q = ±37r/4, while the 
contribution of Y^{uJi\ q,u) relative to that of F' 1 '^; q,co) increases in the a polarization 
due to the polarization factor. 

At q = 37r/4, the contribution of Y^ 2 \u^;q,uj) becomes even larger than that of 
yd) 

(u^;q,u). Since the contribution of FW(w-;g,w) in the a polarization could be esti- 
mated experimentally from the spectra in the ir polarization by multiplying the polarization 
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FIG. 4: Schematic view of the scattering geometry. The scattering plane contains the b and c axes. 
The angle between the incident and scattered x-rays is fixed 130 degrees. 

factor, we may confirm from experiment the large contribution of Y^(u^',q,u) in the a 
polarization at q = 37r/4. We also expect to observe the difference in the spectral shapes 
between in the tt and a polarizations, since the spectral shape of Y^ 2 \uj^] q, uS) has more 
weights at higher to than that of Y^(uf;q,u). Within the present finite-size calculation, 
the spectral peak shifts to a higher energy position in the a polarization from the position 
in the n polarization. Although the system of 16 spins may be too small to discuss spec- 
tral shapes in detail, the fact that more weights exist at higher to in Y^(tOi;q,uj) than in 
y(i) 

(cji, q,u) has been known from the more precise calculation based on the Bethe Ansatz 
solution,— if f^\oji) and ff\uJi) are neglected (actually they are negligible). 

IV. CONCLUDING REMARKS 

We have studied the magnetic excitations in the L-edge RIXS in one-dimensional un- 
doped cuprates. We have analyzed the second-order dipole allowed process through the 
intermediate state, in which there is no spin degree of freedom at the core-hole site. This 
nature of the intermediate state is found to affect strongly the transition amplitudes of spin 
excitations not only at the core-hole site but also at neighboring sites in the final state. 
This tendency is found to increase with decreasing values of Y. The spherical symmetry in 
the spin space in the ground state makes our analysis transparent, making it possible to 
analyze the amplitudes giving rise to excitations not only on the core-hole site but also on 
the neighboring sites. We have evaluated such amplitudes in a finite-size ring. Note that 
the analysis of the RIXS from the antiferromagnetic ordered state was complicated due to 
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q=3u/4 q=-3:r/4 

— , — , — i — , — , — , — | 0.03 | — , — , — , — i — , — , — ,— 




-3-2-1 -3 -2 -1 

to/J (o/J 



FIG. 5: (Color online) RIXS spectra as a function of energy loss u for the momentum transfer 
projected on the b axis q = 37r/4 and — 3tt/4. Black thick solid, red thick broken, and blue thin 
solid lines are total spectrum, Y^(u>i;q,uj), and (ujf, q, oj), respectively. The is set to give 
rise to the peak in the absorption spectra. J = 260 meV and T/J = 1.2. The calculated spectra 
are convoluted with the Lorentzian function with the half width of half maximum 78 meV. 

the presence of the asymmetry in the spin space.— The RIXS spectra have been expressed 
as the one-spin correlation function in the channel with changing polarization and as the 
two-spin correlation function in the channel without changing polarization. We have demon- 
strated that the contribution of the two-spin correlation function could be observed from 
the polarization analysis to possible RIXS spectra in Sr 2 Cu0 3 . 
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FIG. 6: Absorption coefficient A{ui) as a function of photon energy cjj. T/J = 1.2. The origin of 
energy is set to correspond to Ui = e c 



r -corc • 



Appendix A: Absorption coefficient 

The L 2! 3-absorption coefficient Aj(ui) (j = 3/2 or j = 1/2) close to 2p — > 3d x 2_ y 2 
transition may be given by the formula, 

M^) «EiWo)] 2 fM , f f r/7r fM2 , r2 - ( A1 ) 

" 1^ + e 9 - e core - e^j + l 

By substituting the eigenvalues and the eigenstates evaluated on finite-size chain into 
Eq. (lAlj) . we obtain Aj(u>i). Figure [6] shows the calculated Aj(u>i) as a function of pho- 
ton energy. The origin of photon energy is set to be Ui = e core , and T/J = 1.2. The 
calculated curve is found very close to the Lorentzian shape. The peak position is slightly 
shifted from uji = e core ; uji = u;° = e core + 0.8 J for Y / J = 1.2. 
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